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Abstract: Climate change may affect water resources by altering various processes in natural ecosystems. 
Dynamic and statistical downscaling methods are commonly used to assess the impacts of climate change 
on water resources. Objectively, both methods have their own advantages and disadvantages. In the 
present study, we assessed the impacts of climate change on water resources during the future periods 
(2020-2029 and 2040-2049) in the upper reaches of the Kaidu River Basin, Xinjiang, China, and 
discussed the uncertainties in the research processes by integrating dynamic and statistical downscaling 
methods (regional climate models (RCMs) and general citculation modes (GCMs)) and utilizing these 
outputs. The reference period for this study is 1990—1999. The climate change trend is represented by 
three bias-corrected RCMs (i.e., Hadley Centre Global Environmental Model version 3 regional climate 
model (HadGEM3-RA), Regional Climate Model version 4 (RegCM4), and Seoul National University 
Meso-scale Model version 5 (SUN-MM5)) and an ensemble of GCMs on the basis of delta change 
method under two future scenarios (RCP4.5 and RCP8.5). We applied the hydrological SWAT (Soil and 
Water Assessment Tool) model which uses the RCMs/GCMs outputs as input to analyze the impacts of 
climate change on the stream flow and peak flow of the upper reaches of the Kaidu River Basin. The 
simulation of climate factors under future scenarios indicates that both temperature and precipitation in 
the study area will increase in the future compared with the reference period, with the largest increase of 
annual mean temperature and largest percentage increase of mean annual precipitation being of 2.4°C and 
38.4%, respectively. Based on the results from bias correction of climate model outputs, we conclude that 
the accuracy of RCM (regional climate model) simulation is much better for temperature than for 
precipitation. The percentage increase in precipitation simulated by the three RCMs is generally higher 
than that simulated by the ensemble of GCMs. As for the changes in seasonal precipitation, RCMs exhibit 
a large percentage increase in seasonal precipitation in the wet season, while the ensemble of GCMs 
shows a large percentage increase in the dry season. Most of the hydrological simulations indicate that the 
total stream flow will decrease in the future due to the increase of evaporation, and the maximum 
percentage decrease can reach up to 22.3%. The possibility of peak flow increasing in the future is 
expected to higher than 99%. These results indicate that less water is likely to be available in the upper 
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reaches of the Kaidu River Basin in the future, and that the temporal distribution of flow may become 
more concentrated. 
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1 Introduction 


Global climate change may affect water resources worldwide by altering various processes in 
natural ecosystems (Xu and Singh, 2004; Mupenzi et al., 2011). Landman (2007) estimated that 
greenhouse gas emissions are expected to cause global warming to continue throughout the 21‘ 
century. However, the impacts of climate change will vary between regions because of variations 
in the circulation of the atmosphere and the related moisture (Lana et al., 2016). Therefore, having 
reliable information related to climate change at regional scales is crucial to providing effective 
adaption and mitigation measures to climate change. Currently, studies on the impacts of climate 
change primarily concentrates on using general circulation models (GCMs). The coarse resolution 
of GCMs makes it difficult to provide precise information at a regional scale. Statistical 
downscaling method (SDM) and dynamic downscaling method (DDM) are two approaches that 
have been developed for downscaling the climate variables from GCMs (Mearns et al., 1999; 
Busuioc et al., 2006). 

The SDM and DDM use different techniques to extract climate change signals. The SDM 
predicts future regional climate by establishing a statistical relationship between GCM (general 
circulation model) outputs and multi-year observational climate data at large and regional scales. 
However, this approach does not produce detailed spatial information (Christensen and 
Christensen, 2007; Zhang and Huang, 2013). In contrast, the DDM can produce climate 
projections with higher spatial resolution using a regional climate model (RCM) coupled to GCM. 
This approach is heavily influenced by the boundary conditions provided by GCMs. The SDM 
has the primary advantages of being computationally inexpensive and easy to apply, whereas 
running an RCM requires a large amount of computational resources and extensive knowledge of 
computing. Therefore, SDMs are commonly considered as complementary to DDMs (Sachindra 
et al., 2014). 

To date, several RCM-coordinated projects have been developed, including Prediction of 
Regional Scenarios and Uncertainties for Defining European Climate Change Risks and Effects 
(PRUDENCE; Christensen et al., 2007), ENSEMBLES (Hewitt, 2004), Regional Climate Model 
Intercomparison Project (RMIP; Fu et al., 2005), and Coordinated Regional Downscaling 
Experiment (CORDEX; Giorgi et al., 2009), etc. However, most of the domain covered by 
RCM-coordinated projects is European and American continents, and only a few outputs are 
available in East Asia (Coordinated Regional Downscaling Experiment-East Asia (CORDEX-EA) 
being one of them). The CORDEX project (http://www.cordex.org/), sponsored by the World 
Climate Research Programme, aims to overcome the deficiencies of Coupled Model 
Inter-comparison Project Phase 5 (CMIP5) GCMs by using different RCMs at a high resolution 
over twelve different domains (Giorgi et al., 2009). The domain of CORDEX-EA covers East 
Asia, India, western Pacific Ocean, and northern Australia. The number of grid points is 233 
(west-east) by 197 (north-south) horizontal grids centering at 22°02’N and 118°58’E with a 
nominal horizontal resolution of 50 km. The Met Office Hadley Centre developed the Hadley 
Centre Global Environmental Model, version 2-Atmosphere-Ocean (HadGEM2-AO) GCM, and 
the GCM was provided by the National Institute of Meteorological Research, Korea 
Meteorological Administration. All RCMs used the HadGEM2-AO GCM as boundary conditions 
under two scenarios (RCP4.5 and RCP8.5) from the Intergovernmental Panel on Climate Change 
(Bellouin et al., 2011). 

In arid and semi-arid mountainous areas with fragile ecological conditions, water resources 
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mainly come from glaciers, snowmelt and precipitation, and are vulnerable to the negative 
impacts of climate change. The Tarim River is situated in the south of the Tianshan Mountains 
and is the longest inland river in China. The Kaidu River, located in Bayingol Mongolian 
Autonomous Prefecture, Xinjiang of China, serves as a main tributary of the Tarim River. The 
source region of the Kaidu River connects the Bayanbulak Swan Nature Reserve and eventually 
flows into the Bosten Lake, the largest inland freshwater lake in China. The surrounding alpine 
mountains provide abundant water sources to the grassland and irrigation land in the downstream 
region. Therefore, the Kaidu River is important for maintaining the stability of the environment as 
well as the productivity and livelihood of the local residents. Due to arid and semi-arid 
environment, the ecosystem in the Kaidu River Basin is especially fragile and sensitive to climate 
change and human activities. Thus, it is necessary to determine the hydrological responses of the 
Kaidu River Basin within the context of climate change. Numerous studies have been conducted 
to analyze the impacts of climate change on water resources in the Kaidu River Basin (Yang and 
Cui, 2005; Xu et al., 2008; Jiang et al., 2011; Liu et al., 2011; Wang et al., 2013). These studies 
utilized the outputs of GCMs coupled with hydrological models, and analyzed the hydrological 
responses of the river basin based on the predicted climate change. For example, Liu et al. (2011) 
analyzed the outputs from 36 GCMs in two hydrological models, and revealed that the snowmelt 
rate increases and the snowmelt period will be pushed forward in time via changes in temperature. 
Fang et al. (2015) compared the application of five precipitation and three temperature correction 
methods for downscaling RCM simulations over the Kaidu River Basin, and pointed out that 
different correction methods for precipitation could create large differences in the results, while 
correction methods for temperature showed relatively consistent results. Xu et al. (2016) used two 
SDMs to downscale the Hadley Centre Coupled Model version 3 (HadCM3) outputs in order to 
project the future climate scenarios in the Kaidu River Basin. Under scenarios A2 and B2, they 
found that the basin would experience an increasing trend in temperature as well as an indistinct 
change trend in precipitation during 2011—2040. 

From the above studies, we can see that the impacts of climate change on water resources in the 
Kaidu River Basin at different temporal and spatial scales remain an open research question. 
Selecting appropriate and optimal climate and hydrological models is the key to solving this 
problem. Conducting a comprehensive assessment of hydrological responses to climate change 
using statistical or dynamical downscaling GCM outputs in conjunction with hydrological models 
has been attempted in some previous studies (Wilby et al., 2000; Fu et al., 2011). In the present 
study, outputs from different dynamic downscaling methods under the same boundary conditions 
were used to assess the impacts of climate change on water resources in the upper reaches of the 
Kaidu River Basin. It should be noted that we used the multi-GCM set averaging approach as a 
complement to avoid the systematic bias that originates from the single boundary condition. 

Reducing the uncertainties associated with the downscaling methods and determining the 
structure of hydrological models are the major issues that to be addressed to accurately estimate 
the impacts of climate change on water resources. DDMs can produce high-resolution climate 
projections by nesting a GCM and an RCM. However, previous applications of RCMs have 
mainly obtained initial boundary conditions from a single GCM. It is possible that the use of a 
single GCM amplifies the uncertainty. Moreover, a large number of the uncertainties introduced 
in RCMs are caused by the internal dynamics of RCMs as well as by external forcing. Thus, 
applying the outputs from a single RCM in a hydrological model may potentially amplify the 
uncertainty (Foley, 2010; Wehner, 2013). With regard to the structure of hydrological models, the 
distributed hydrological models have more advantages when compared with lumped conceptual 
models. First, a distributed hydrological model can reflect the physical mechanisms of 
hydrological processes. This produces less uncertainty than a lumped conceptual model. Second, 
a distributed hydrological model that has been validated is applicable to areas lacking 
meteorological data. 

For this study, one of our goals is to reduce the uncertainty associated with the use of a single 
downscaling method in hydrological modeling. Therefore, we used three RCMs (i.e., Hadley 
Centre Global Environmental Model version 3 regional climate model (HadGEM3-RA), Regional 
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Climate Model version 4 (RegCM&4), and Seoul National University Meso-scale Model version 5 
(SUN-MMS)) and an ensemble of GCMs under two scenarios (RCP4.5 and RCP8.5) to project 
the trends of climate change in the upper reaches of the Kaidu River Basin. A distributed 
hydrological model was then used to analyze the hydrological responses of the Kaidu River Basin 
to climate change during future period (2020-2049). 


2 Materials and methods 


2.1 Study area 


The Kaidu River Basin is located in Xinjiang, China and lies along the central and southern slopes 
of the Tianshan Mountains in Xinjiang (42°14'N—-43°21'N, 82°58'E-86°05'E; Fig. 1). The basin 
covers an area of 18.9x10° km? above the Dashankou hydrological station (namely the upper 
reaches of the Kaidu River Basin). The Kaidu River has a total length of 560 km. The origin of the 
river is in Yilianhabierga Mountain. The river flows through the Yulduz and Yanqi basins, and 
eventually enters into the Bosten Lake. The Kaidu River at the Dashankou hydrological station 
contributes 35.1x108 m? in annual runoff volume or 85% of the annual inflow into the Bosten Lake 
(Zhang et al., 2014). The elevation of the upper reaches of the Kaidu River Basin ranges from 1293 
to 4796 m a.s.l. (Fig. 1), with an average elevation of 3100 m a.s.l. (Liu et al., 2017). 

In this study, one meteorological station (Bayanbulak) and two hydrological stations 
(Bayanbulak and Dashankou) were utilized to document climate conditions in the upper reaches 
of the Kaidu River Basin. Data from another meteorological station (i.e., Baluntai) near this basin 
were also employed in this research (Fig. 1). The upper reaches of the Kaidu River Basin features 
a temperate continental arid climate, with mean annual precipitation of approximately 270 mm. 
The study area has a long winter that lasts from November to March next year. Annual mean 
temperature and the extreme minimum temperatures are —3.85°C and —48.1°C, respectively (Xue 
et al., 2016). 
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Fig. 1 Overview and the location of the upper reaches of the Kaidu River Basin 


2.2 Climate model outputs and bias correction 


In this study, both dynamic and statistical downscaling outputs were used to explore the impacts 
of climate change on water resources. The outputs of three dynamic downscaling models were 
obtained from RCMs that were used in the CORDEX-EA project. In order to obtain more 
representative statistical downscaling outputs, we used the delta change method to calculate the 
relative change in variables (Onyutha et al., 2016). First, we examined the outputs of all GCMs 
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from Coupled Model Inter-comparison Project Phase 5 (CMIP5) by comparing the GCM 
historical data with observed records in the study area. Moreover, we then excluded the simulated 
temperature and precipitation data with the standard deviations (between GCM runs and observed 
data) being larger than 1. Eventually, 30 GCMs were selected and the relative change in variables 
was calculated. The daily temperatures for future scenarios were obtained by adding the observed 
daily series (1990-1999) to the differences in mean monthly temperature between GCM scenario 
series and reference series. The simulations for daily precipitation were generated via multiplying 
the observed daily series (1990-1999) by the rate of increase in the mean monthly precipitation in 
two future periods (2020-2029 and 2040-2049). A summary of RCMs and GCMs used in the 
present study is given in Table 1. 


Table 1 Regional climate models (RCMs) and general circulation models (GCMs) used in the present study 


RCM/GCM Resolution GCM Scenario Variable Time period Reference 
HadGEM3-RA 0.44° HadGEM2-AO tee Pre, Tmax, Tmin ER E Davies et al. (2005) 
RegCM4 50km HadGEM2-AO ee Pre, Tmax, Tmin Ra ph setae Giorgi et al. (2012) 
SUN-MMS 50km  HadGEM2-AO ee Pre, Tmax, Tmin snes Ee Cha and Lee (2014) 
GCMs 1.4°—5.0° CMIP5 been Pre, Tmax, Tmin Ao etnias, Taylorer et al. (2011) 


Note: RCM, Regional climate model; GCM, general circulation model; HadGEM3-RA, Hadley Centre Global Environmental Model 
version 3 regional climate model; RegCM4, Regional Climate Model version 4; SUN-MMS, Seoul National University Meso-scale 
Model version 5; HadGEM2-AO, Hadley Centre Global Environmental Model, version 2-Atmosphere-Ocean; Pre, precipitation; Tmax, 
daily maximum near-surface air temperature; Tmin, daily minimum near-surface air temperature. 


When using the delta change method, we removed the GCM with a large bias in temperature 
and precipitation simulation; therefore, we did not consider the bias in SDMs in this study. 
Directly utilizing the outputs of RCMs in a hydrological model is not recommended because the 
raw RCM data have systematic errors (Wilby et al., 2000). Therefore, we used two methods to 
correct the outputs of RCMs. 

The local intensity scaling (LOCI) method was used for precipitation correction, and the linear 
scaling (LS) method was applied for temperature correction. The correction of precipitation data 
by using LOCI method involves three steps (Schmidli et al., 2006). First, we calculated a wet-day 
threshold from the raw precipitation series by ensuring that the threshold exceedance matches the 
wet-day frequency of observed precipitation series. Second, a scaling factor was calculated and 
then used to ensure that the mean of simulated precipitation is equal to that of observed 
precipitation. Third, we employed these thresholds and scaling factors based on the reference 
period (1990-1999) to adjust precipitation for the future (2020—2029 and 2040-2049) predictions. 

The LS method assumes that biases are time-invariant (Lenderink et al., 2007). Moreover, the 
method also assumes that the monthly mean of raw temperature data is the same as that of 
observed data. We adjusted temperature data via monthly correction values based on the 
differences between observed and raw data for reference (1990-1999) and future (2020-2029 and 
2040-2049) periods. These adjustments were done with the following equations: 

Peon = Te I), (1) 
Reo. =f + (Tp -T7 ), (2) 
where, T(cor);" and 7(cor);" are the corrected temperatures on the i‘ day of the reference period 
and future period (°C), respectively; T;" and T;" are the simulated temperatures on the i" day of 


the reference period and future period (°C), respectively; and T?’ and T represent the 


mean value of observed and simulated temperatures (°C), respectively, during the reference 
period at a given month m. 

After bias correction, the outputs from RCMs can be used to analyze the impacts of climate 
change on water resources and extreme events. The present study did not apply the outputs of 
SDMs to extreme events. This is because the delta change method only accounts for change in the 
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mean while ignores other possible types of change in the distribution of the variables. 

We analyzed the future (2020-2029 and 2040-2049) change of peak flow in the upper reaches 
of the Kaidu River Basin by comparing reference period (1990-1999) to the outputs of three 
RCMs that were used to simulate future stream flow. The P-III frequency curve (Geng, 2002) and 
percentile value method (Dimitropoulos et al., 2009) were used to produce the frequency curve 
for peak flow of the upper reaches of the Kaidu River Basin under RCP4.5 and RCP8.5 scenarios 
during the periods 2020-2029 and 2040-2049. The modeling time of SUN-MM5 under RCP8.5 
scenario was from 2020 to 2035. Therefore, only the period 2020-2029 was considered in the 
frequency curve of peak flow in the upper reaches of the Kaidu River Basin under RCP8.5 
scenario. 


2.3 Hydrological assessment 


We assessed the impacts of climate change on water resources in the upper reaches of the Kaidu 
River Basin by using outputs from three RCMs and an ensemble of GCMs to drive a hydrological 
model. Instead of directly applying the outputs of RCMs, we used the historical simulation data 
and the observational data to find the bias scaling factors, and then corrected the bias of future 
outputs from RCMs accordingly. In addition, we extracted the climate change signals of GCM 
data based on its historical data and future simulations. After input data processing and 
hydrological model constructing, the impacts of climate change on water resources were 
investigated during a reference period (1990-1999) and two future periods (2020-2029 and 
2040-2049) under RCP4.5 and RCP8.5 scenarios in the upper reaches of the Kaidu River Basin. 
Detailed steps are shown in Figure 2. 
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Fig.2 Framework of the research methodology for hydrological assessment of the impacts of climate change on 
water resources in the upper reaches of the Kaidu River Basin. RCM, Regional climate model; GCM, general 
circulation model; SWAT, Soil and Water Assessment Tool. 


A physically based, semi-distributed, and temporally continuous watershed model called the 
Soil and Water Assessment Tool (SWAT) was used for stream flow simulation in the present study. 
The key processes of hydrological simulation in SWAT consist of precipitation, vegetation 
interception, discharge, groundwater, snowmelt, actual/potential evapotranspiration and so on. 
More details on the hydrological processes and model information can be found in SWAT 
manuals (http://www.brc.tamus.edu/). The processes of snowmelt and evapotranspiration are 


BA Wulong et al.: Simulating hydrological responses to climate change using dynamic and statistical... 


especially important to the local hydrological cycle and dominate the spatial and temporal 
variations of water resources. 

In this study, the SWAT model was established using topographical, land use, soil map and 
weather data. A set of 90 mx90 m topographical maps was obtained from the Shuttle Radar 
Topography Mission (SRTM) of the United States Geological Survey. Land use data were 
obtained from the Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences. 
Soil map data at a scale of 1:1,000,000 were derived from the Harmonized World Soil Database. 
Daily weather data, including precipitation and maximum/minimum temperature at Bayanbulak 
and Baluntai meteorological stations, were obtained from the China Meteorological Data Sharing 
Service System (http://cdc.cma.gov.cn/). Daily streamflow dataset from the Dashankou 
hydrological station was provided by the Xinjiang Tarim River Basin Management Bureau. 


2.4 SWAT model performance evaluation 


We employed Nash and Sutcliffe efficiency (NS), percent bias (PBIAS) and correlation 
coefficient (R?) as hydrological model evaluation statistics (Eqs. 3-5, respectively; Moriasi et al., 
2007; Rahman et al., 2014): 


NS=1- 2 aaa (3) 

> Qn- 

a 2 

B Qs x 0,4) 
PBIAS(%) =| == — x 100%, A 
0; 
T — —_— 2 
Yor -2,)Q,, -2,) 

(5) 


£o, -0 |e. -a f” | 


where, Qo; and Qs are the observed and simulated variables at time t, and Q, is the mean of 


observed variables. 

NS indicates the strength of the relationship between observation and simulation, with the 
values ranging from —œ to 1.0. The higher the value (closer to 1.0), the better the model 
performance. PBIAS indicates the average tendency of simulated data compared with observed 
data. Positive and negative values indicate an overestimation and an underestimation of 
observation, respectively. R? represents the proportion of total variance in the observed data that 
can be explained by the simulated data. The performance of model is considered to be improved 
as R? approaches to 1.0. 


3 Results and discussion 


3.1 Hydrological modeling performance 


Hydrological modeling provides critical information for the assessment of water resources 
(Sivapalan et al., 2003). The use of unreasonable parameter sets might modify model responses in 
the future scenarios. Our SWAT model was calibrated and validated via the daily observed flow 
data at the Dashankou hydrological station. Sequential Uncertainty Fitting (SUFI-2) was used for 
model parameter sensitivity analysis. Twenty-seven hydrological parameters were evaluated, but 
only ten parameters were selected to avoid the problems related to over-parameterization. The 
three most sensitive parameters were the base flow alpha factor, the precipitation lapse rate, and 
the temperature lapse rate. This means that the base flow, elevation band of precipitation and 
temperature had the greatest impact on stream flow. The parameters were then manually 
calibrated. The base flow alpha factor was adjusted to 0.78 from the default value of 1.00. The 
temperature and precipitation lapse rates were adjusted to —1.45°C/km and 137 mm/km, 
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respectively. Figure 3 shows the performance of SWAT model during calibration (1986-1995) and 
validation (1996-2005) periods, and a summary of the results is shown in Table 2. The results of 
NS, PBIAS, and R? indicate that the model exhibits a good performance during calibration and 
validation periods. 
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Fig. 3 Calibration (a) and validation (b) results for daily stream flow during periods 1986-1995 and 1996-2005 
at the Dashankou hydrological station, respectively 


Table 2 Model performance rating based on Moriasi et al. (2007) and Rahman et al. (2014) 


Period NS PBIAS (%) R 
Calibration (1986-1995) 0.69 1.20 0.72 
Validation (1996-2005) 0.65 1.43 0.67 


3.2 Bias correction of climate model 


In this study, three RCM-generated variables were corrected via bias correction during the 
reference period (1990—1999) using observed data from two meteorological stations (Baluntai and 
Bayanbulak). The LOCI and LS bias correction methods were applied to the output variables of 
temperature and precipitation. After correcting the historical period, the same correction factors 
were applied to the RCP4.5 and RCP8.5 future scenarios. We evaluated the correction 
performance by generating Taylor diagrams that compared the raw and corrected RCM variables 
with the observations. The performances of bias corrections compared with observations are 
presented in terms of R?, standard deviation (SD), and RMSE (root mean square error). The 
Taylor diagram has three axes: R? (azimuthal position), SD (x and y axes) and centered RMSE 
(Fig. 4). The observed value as a reference (expressed as REF) intersects the x axis. 

Originally, the R? values of raw precipitation ranged from 0.14 to 0.94. RegCM4 and 
SUN-MM5S exhibited the highest and lowest coefficient values, respectively. The R? values of the 
maximum and minimum temperatures were similar, with the values ranging from 0.94 to 0.97. 
RegCM64 presented the best SD performance for raw precipitation, while HadGEM3-RA had the 
largest SD value. From the results shown in Figure 4, we can infer that the accuracy of RCM 
simulation is much better for temperature than for precipitation. The bias correction results during 
the reference period (1990-1999) for RCM outputs were satisfactory. All coefficient values of 
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corrected variables were higher than 0.99. The R? and RMSE values of corrected RCMs were also 
very close to those of field observations. These results suggest that the effect of bias correction is 
efficient and acceptable. 
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Fig. 4 Taylor diagrams for the raw regional climate model (RCM) results of monthly precipitation (a), minimum 
temperature (b), and maximum temperature (c) versus the corresponding bias corrected results. The Taylor diagram 
has three axes: correlation coefficient (R?; azimuthal position, black dotted line), standard deviation (SD; radical 
distance, x and y axes) and root mean square error (RMSE; distance from REF (reference) point, yellow line). The 
REF on the x axis represents observation. HadGEM3-RA, Hadley Centre Global Environmental Model version 3 
regional climate model; RegCM4, Regional Climate Model version 4; SUN-MM5S, Seoul National University 
Meso-scale Model version 5. A, B and C represent the results of HadGEM3-RA, RegCM4 and SUN-MMS models, 
respectively. D, E and F represent the bias corrected results of HadGEM3-RA, RegCM4 and SUN-MMS5 models, 
respectively. 


3.3 Future trends of climate change 
3.3.1 Changes in seasonal temperature and precipitation 


To investigate climate change in different seasons, we compared the changes of seasonal 
temperature and precipitation (percentage change) from the ensemble of GCMs and three RCMs 
fe during the reference period (1990-1999) and future periods (2020-2029 and 2040-2049) under 
ph RCP4.5 and RCP8.5 scenarios (Fig. 5). In the present study, we defined the wet (from April to 
i= September) and dry (from October to March) seasons to analyze the impacts of seasonal climate 
change. 

All the outputs from the three RCMs and the ensemble of GCMs show that the temperature in 
the studied region (i.e., the upper reaches of the Kaidu River Basin) will continue to rise during 
2020-2049 under both scenarios (RCP4.5 and RCP8.5). The increase in temperature will be more 
prominent under RCP8.5 scenario than under RCP4.5 scenario. The results confirm that the 
greenhouse gas emission level directly affects the temperature change (Solomon et al., 2007). 
Figure 5 clearly depicts that the incremental changes in the maximum and minimum temperatures 
of three RCMs and the ensemble of GCMs are similar within each individual season. However, 
the variation in the maximum and minimum temperatures varies between seasons for different 
models. The maximum temperature changes of RegCM4, SUN-MM5 and ensemble of GCMs will 
increase more significantly in the dry season than in the wet season, excluding values from 
HadGEM3-RA. The minimum temperature changes in RegCM4 and SUN-MM5 will also 
increase more dramatically in the dry season then in the wet season. 

The simulation of precipitation by the three RCMs and the ensemble of GCMs shows that 
percentage changes of precipitation in the dry and wet seasons are obviously increasing under the 
two climate change scenarios (RCP4.5 and RCP8.5), with the exceptions of SUN-MM5 under 
RCP4.5 and HadGEM3-RA under RCP8.5 (Fig. 5). It should be noted that the precipitation 
increases modeled by the three RCMs are more pronounced in the wet season than in the dry 
season. In contrast with RCMs, the increments of precipitation from the ensemble of GCMs are 
smaller in the wet season than in the dry season. 
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Fig. 5 Changes in temperature (a-d) and percentage changes in precipitation (e-f) in the dry and wet seasons 
from the ensemble of GCMs and three RCMs during the future periods (2020-2029 and 2040-2049) under 
RCP4.5 and RCP8.5 scenarios compared with the reference period (1990-1999) 


3.3.2 Changes in annual temperature and precipitation 


Relative to the reference period (1990-1999), the simulated outputs from all models show a 
continuous increase in annual mean temperature for the future periods (2020-2029 and 
2040-2049) under both RCP4.5 and RCP8.5 scenarios (Fig. 6a). The simulated changes in 
magnitude of annual mean temperature under RCP8.5 scenario reach up to 2.4°C by the mid-21* 
century. The simulated percentage changes in mean annual precipitation also show an upward 
trend during the future periods compared with the reference period (Fig. 6b). The result is 
consistent with Chen et al. (2014) that the climate in northwestern China is gradually changing 
from warm dry to warm wet in the future. 

In contrast to the results of annual mean temperature simulations, the mean annual precipitation 
values simulated by the three RCMs and the ensemble of GCMs exhibit obvious variations. The 
percentage increase in mean annual precipitation from the three RCMs is larger than that from the 
ensemble of GCMs. For example, the simulated percentage change in precipitation from RegCM4 
reaches up to 38.4% under RCP8.5 scenario, while the value is only 5.0% from the ensemble of 
GCMs. 
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Fig. 6 Changes in annual mean temperature (a) and percentage changes in mean annual precipitation (b) from 
the ensemble of GCMs and three RCMs under RCP4.5 and RCP8.5 scenarios during the future periods 
(2020-2029 and 2040-2049) 


3.3.3 Inter-annual variation trends of extreme precipitation 


In the context of human-caused global warming, water circulation is predicted to accelerate 
(Handel and Risbey, 1992). Extreme precipitation, floods, and other natural disasters have serious 
impacts on regional economic development and the lives of local residents. In this study, we 
compared the extreme precipitation data calculated by the percentile value method 
(Dimitropoulos et al., 2009) for three RCMs under RCP4.5 and RCP8.5 scenarios during the 
future period 2020-2049 (Fig. 7). The threshold percentile value for extreme precipitation was 
defined as 95%. Figure 7 illustrates that under both RCP4.5 and RCP 8.5 scenarios, the 
inter-annual variation trends for all simulated extreme precipitation data from the three RCMs 
show a slight upward trend, with the exception of SUN-MM5 under RCP8.5 scenario. This means 
that the risk of flood disasters in the upper reaches of the Kaidu River Basin is likely to rise in the 
next 30 years. 
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Fig. 7 Interannual variation trends of extreme precipitation from the three RCMs under RCP4.5 and RCP8.5 
scenarios during the future period 2020-2049 


3.4 Hydrological responses to climate change 


3.4.1 Changes in annual stream flow under different scenarios 


Figure 8 shows the percentage changes of annual stream flow in the upper reaches of the Kaidu 
River Basin as simulated by hydrological models according to bias corrected outputs of three 
RCMs and the ensemble of GCMs under RCP4.5 and RCP8.5 scenarios during the two future 
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periods (2020-2029 and 2040-2049). The results of simulations from the three RCMs and 
ensemble of GCMs indicate that the annual stream flow will tend to decrease during 2020-2049 
under the two scenarios, with the exceptions of RegCM4 under RCP8.5 scenario and SUN-MM5 
under both scenarios. The simulations from HadGEM3 and the ensemble of GCMs show a higher 
annual stream flow under RCP4.5 scenario than under RCP8.5, while RegCM4 and SUN-MM5 
simulate less annual stream flow under RCP4.5 scenario than under RCP8.5. As shown in Figure 
8, the largest percentage decrease in annual stream flow is 20.1% and 22.3% during the period 
2040-2049 under RCP4.5 and RCP8.5 scenarios, respectively, which are simulated by the 
ensemble of GCMs. From the abovementioned results, we can conclude that the changes in 
temperature and precipitation could lead to a remarkable change in snowmelt and evaporation, 
ultimately affecting stream flow in the upper reaches of the Kaidu River Basin. First, as the 
temperature rises, the amount of evaporation will inevitably rise, causing a large change in runoff. 
Second, in the long run, glacial reserves will decrease as the temperature rises. After the river of 
glacial meltwater reaches the maximum stream flow, there will be a gradual declining trend. 
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Fig. 8 Percentage changes in annual stream flow in the upper reaches of the Kaidu River Basin from the 
ensemble of GCMs and three RCMs under RCP4.5 and RCP8.5 scenarios during the future periods (2020-2029 
and 2040-2049) 


3.4.2 Frequency curves of peak flow in the upper reaches of the Kaidu River Basin 


During the future periods (2020-2029 and 2040-2049), among the 11 frequency curves of peak 
flow in the upper reaches of the Kaidu River Basin that were plotted from the three RCMs under 
RCP4.5 and RCP8.5 scenarios, only one curve in HadGEM3-RA under RCP8.5 scenario during 
the period 2040-2049 suggests a potential decrease in peak flow compared with the reference 
period (Fig. 9). Eight frequency curves plotted by HadGEM3-RA, SUN-MM5, and RegCM4 
simulations indicate that the peak flow has a greater than 99% probability of increasing in the 
future. The other two frequency curves plotted by RegCM4 simulations show that the peak flow 
has 70% and 85% probability of increasing under RCP4.5 scenario during 2020-2029 and 
2040-2049, respectively. 

The frequency curves of peak flow indicate that the peak flow will not decrease despite the fact 
that the stream flow in the upper reaches of the Kaidu River Basin is expected to decline. In other 
words, the temporal distribution of flow becomes more concentrated. The analysis shows that the 
continuous rise in temperature has exacerbated the occurrence of extreme hydrological events in 
the study area, leading to the redistribution of water resources. 
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Fig. 9 Frequency curves of peak flow in the upper reaches of the Kaidu River Basin from the three RCMs 
during the reference period (1990-1999) and the future periods (2020-2029 and 2040-2049) under RCP4.5 and 
RCP8.5 scenarios 


4 Uncertainties analysis 


It is obvious that many uncertainties exist in the assessment of hydrological responses to climate 
change. They mainly include: (1) uncertainties to different GCMs; (2) uncertainties to different 
downscaling methods; (3) uncertainties about different future scenarios; (4) uncertainties to 
different bias correction methods; and (5) uncertainties of different hydrological models. 
Generally speaking, GCM is the largest source of uncertainties in the simulation of future 
hydrological predictions (Chen et al., 2011). In the study of the impact of hydrological forecast 
uncertainties on decision-making, Buytaert et al. (2010) pointed out that the uncertainties of GCM 
accounted for 50.8% of the total uncertainty. In this study, the three RCMs (i.e., HadGEM3-RA, 
SUN-MMS, and RegCM4) projected large differences in temperature and precipitation in the 
future. The increases in temperature during the future periods obtained from the three RCMs vary 
from 0.6°C to 2.0°C under RCP4.5 scenario and from 1.1°C to 2.4°C under RCP8.5 scenario. The 
simulated precipitation also has relative large differences among the three RCMs. Specifically, the 
simulated percentage changes of annual precipitation in the future range from 4.0% to 29.5% 
under RCP4.5 scenario and from 4.3% to 38.4% under RCP8.5 scenario by using the three RCMs. 
We found that the downscaled outputs in different physical parametric methods are distinct even 
at the same boundary condition. Moreover, the raw RCM data usually deviate from the observed 
data. Thus, we conclude that it is necessary to correct the data before hydrological modeling, 
which confirms the results of Ouyang et al. (2014, 2015) regarding the uncertainties to different 
downscaling methods. In this study, we used more confident GCM ensembles as a benchmark to 
avoid the systematic bias. In addition, although SWAT can generally capture the hydrological 
processes in the study area, only 10 most sensitive parameters were selected for calibrating. The 
neglected parameters might potentially relate to the climate change factors, such as the spatial 
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heterogeneity of temperature/precipitation lapse rate, the temporal heterogeneity of channel 
conductivity and so on. 


5 Conclusions 


The outputs of RCMs and GCMs provide consistent results, and the results indicate that the 
temperature and precipitation will increase in the upper reaches of the Kaidu River Basin during 
the future periods (2020-2029 and 2040-2049) when compared with the reference period 
(1990-1999). However, according to the hydrological simulation results generated by RCMs and 
GCMs, the Kaidu River Basin will likely experience a period of less water availability and may 
suffer extensive flood risks in the middle of this century due to the increase of evaporation in the 
future. In contrast with the RCMs, the percentage increases in precipitation from the ensemble of 
GCMs are more pronounced in the dry season than in the wet season. The probability of 
occurrence of peak flow in the future will increase despite the fact that the stream flow in the 
upper reaches of the Kaidu River Basin is expected to decline, which means that the temporal 
distribution of flow becomes more concentrated. 

Our results suggest that the existence of uncertainties in the assessment of hydrological 
responses to climate change likely relates to the uncertainty of the input data, the uncertainty of 
the future climate scenarios, and the uncertainty in the hydrological models. By properly applying 
the bias correction, the confidence of scenario data can be improved. Since the potential 
uncertainties originate from multiple sources, we would like to extend our future study on 
quantifying the uncertainties based on the entire impact analysis procedures. 
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